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ABSTRACT 

We report interferometric observations of the high-mass star-forming object IRAS 23033+5951. 
Our observations reveal two massive molecular cloud cores, designated IRAS 23033+5951-MMS1 and 
IRAS 23033+5951-MMS2. MMS1 has already formed a massive protostar and MMS2 appears to be 
on the verge of doing so. The latter core may be an example of a massive analogue to a "Class 0" star- 
forming object. The more evolved core shows some evidence of N2H + destruction near the protostar, 
consistent with similar findings in low-mass star-forming objects. In addition to the already-known 
prominent HCO + outflow, our SiO 2-1, and CH3OH 2-1 maps show evidence for two more candidate 
outflows, both presumably less powerful than the main one. Both cores are embedded in an elongated 
feature whose major axis is oriented almost exactly perpendicular to the axis of the most prominent 
outflow in the region. Although it has many of the characteristics of a disk, the 87,000 AU (0.42 pc) 
diameter of this structure suggests that it is more likely to be the flattened, rotating remnant of the 
natal molecular cloud fragment from which the star-forming cores condensed. We conclude that IRAS 
23033+5951 is an excellent example of massive star formation proceeding in relative isolation, perhaps 
by the method of monolithic collapse and disk accretion. 

Subject headings: stars: formation — ISM: individual (IRAS 23033+5951) — ISM: molecules — ISM: 
kinematics and dynamics 



1. INTRODUCTION 

\& \ Two theorie s of massive star format ion dominate the 
CN| ■ literature (see iMcKee fc Ost rikcr 2007 for a recent com- 
' prehensive review). According to the first, massive 
[ stars form by the gravitational collapse of molecular 
. cloud cores, which have masses of tens to hundreds 
y—i ' of solar masses. The molecular cloud cores are them- 
t — I selves produced by the fragmentation of molecular cloud 
. clumps, which have masses of hundreds to thousands 
' of solar masses. In this "monolithic collapse" model, 
• t-H I massive cores are envisioned to collapse to form single 
massive stars or small multiple systems via a process 
schematically similar to that which forms low-mass stars 
. <3 i (IMcKee fc Tanl [20031: iKrumholz. McKee. fc KlehU [20051 : 
lKrumholdl2006 ) . The other prominent theory of massive 
star formation, known as competitive accretion, holds 
that all stars begin as low-mass protostars which then 
compete to accrete some fraction of the mass of their par- 
ent mol ecular cloud clump, sometimes coalescing in the 
proces s (IBonnell. Bate, fc Zinnec kcr 1998t IBonnell et aT] 
I2001allbl : iBonnell. Vine, fcl lBatef2004D . In the compet- 
itive accretion scenario, massive stars are simply those 
which accrete most successfully. It is hard to imagine 
that competitive accretion does not occur at some level 
in all clustered star-forming regions. However, recent 
theoretical work suggests it i s not the dominant mode by 
which massive stars form ([Krumholz. McKee. fc Klein! 
2005). Direct evidence capable of distinguishing these 
modes of star formation remains elusive. To prove that 
massive stars can form by the collapse of individual 
molecular cloud cores would require several conditions to 
be verified by direct observation: (a) a massive protostar 



forming inside a well-defined molecular cloud core, (b) 
an infall signature in that molecular cloud core which is 
consistent with a high rate of protostellar accretion, (c) 
a massive circumstellar disk (whose properties are not 
yet well constrained), and (d) a collimated outflow with 
a high entrained mass. 

Relatively few massive star-forming objects with 
most or all of these characteristics have been identi- 
fied. Perhaps the b e st exa mple is IRAS 20216+4104 
([Cesaroni et al.l 119971 . [1999) , for which there is strong 
evidence of a collapsing Keplerian disk around a massive 
protostar. Other examples, such as Cepheus A HW2 
(jPatel et al.ll2005l ) , are hotly disputed. The principal ob- 
stacle in the study of massive protostellar systems is lack 
of spatial resolution. Because massive protostars typi- 
cally lie several kiloparsecs or more from the Sun, it is 
extremely challenging to spatially resolve their accretion 
disks, even with modern millimeter and submillimeter 
interferometers. Our inadequate understanding of the 
chemistry of massive protostellar systems poses another 
significant challenge, making it difficult to select appro- 
priate molecular transitions with which to study the kine- 
matics of such systems. 

IRAS 23033+5951 (IRAS 23033 hereafter) is an ac- 
tive star-forming region lying at a distance of 3.5 kpc 
in the Cepheus star-forming complex. Evidence of 
ongoing massive star formation has been detected 
in the form of comp a ct radio cont i nuum emission 
(ISridharan et all 120021: iBeuther et alj l2002b[) . maser 
emission (Bcuther et al.1 l2002bD , and two discrete dus t 
continuum sources at 1.2 mm ([Beuther et al.l [2002a). 
However, in contrast to many high-mass star-forming 
regions, relatively few discrete sources have been de- 
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tected a nd there is evidence for only one sub stantial 
outflow (jBeuther et alJl2002d : IMao fe Zengll2003f ). Thus, 
IRAS 23033 provides an attractive opportunity to study 
massive star formation in a relatively quiescent, morpho- 
logically simple environment in which the usual compli- 
cating factors — crowding, abundant contamination from 
radio continuum emission, lack of spatial resolution — are 
mitigated. 

Our goal in studying IRAS 23033 was to characterize 
its kinematics and determine whether it is a candidate 
example of a massive star-forming by monolithic collapse, 
as suggested by its simple morphology and known kine- 
matics. We observed IRAS 23033 in several molecular 
transitions using the Berkeley-Illinois-Maryland Associ- 
ation (BIMA) array. Each spectral line was selected to 
probe a particular kinematic aspect of IRAS 23033 in- 
cluding outflows, rotation, infall, and turbulence. In 
we discuss the observations and the data reduction pro- 
cedure. In Sj3( we discuss the general morphology and 
evolutionary state of IRAS 23033, as determined from 
both the continuum and molecular line observations. In 
$4j we concentrate on the kinematic results derived from 
the spectral line observations. Our results are summa- 
rized in iJS] 

2. OBSERVATIONS AND DATA REDUCTION 

Observations were made over 8 nights between 2003 
October and 2004 May using the BIMA interferometer in 
Hat Creek, CA. All of the observations were made using 
single pointings toward the position a = 23 h 05 m 25 1 ?7, 
5 = +60°08'08" (J2000). Tracks were obtained in the B, 
C, and D configurations of the array and included tunings 
at 86, 89, and 93 GHz. Not all of the spectral lines were 
observed with the same velocity resolution during every 
track. Table [T] provides a complete summary of the lines 
observed, array configurations used, and the spectral line 
velocity resolutions used in most of our analysis. 

The data were flagged, calibrated, CLEANed 
and analyzed using the MIRIAD software package 
(|Sault. Teuben. fe Wrightl 119951 ). Observations of the 
nearby quasars 0102+584 and 2203+317 were made ev- 
ery 30 minutes during the observations and were used 
to correct instrumental and atmospheric phase and am- 
plitude variations in the data. Absolute flux calibra- 
tion was performed using observations of Uranus taken 
during each track. Based on the variations in the flux 
of the quasar from one epoch to the next, we estimate 
the accuracy of our flux calibration to be approximately 
30%. In order to maximize the signal-to-noise ratio of 
the data, we typically used natural weighting of the vis- 
ibilities in generating our images. Where the analysis 
required higher spatial resolution than that provided by 
natural weighting, some of the data were processed us- 
ing a MIRIAD "robust" parameter of 0.5, at the cost of 
a lower signal-to-noise ratio. 

Our 3 mm continuum map was produced by combining 
continuum data from all of the different tunings around 
3 mm. The resulting CLEANed, naturally weighted im- 
age has a beam size of 5'.' 5x4'.' 8 and an rms sensitivity of 
2.9 mJy beam -1 . 

3. MORPHOLOGY: MULTIPLE MASSIVE CORES 

The initial millimeter continuum observations of 
IRAS 23033 with the IRAM 30 m telescope showed a 



singly peaked structure (jBeuther et alj l2002al) . Sub- 
sequent BIMA array observations of the source at 
2.6 m m showed two continuum peaks (jBeuther et al.l 
2002b). Our BIMA array 3 mm continuum image, 
shown in Figure [TJ reveals IRAS 23033 to be composed 
of at least three continuum peaks joined by contigu- 
ous emission, plus a weak fourth peak visible to the 
east at the 3 a level. We designate these three mil- 
limeter continuum peaks as IRAS 23033+595 1-MMS1, 
IRAS 23033+5951-MMS2, and IRAS 23033+5951- 
MMS3 (hereafter, MMS1, MMS2, and MMS3) as shown 
in Figure [TJ We used a simple model to determine the 
fraction of the continuum emission which is associated 
with each of the continuum peaks. The model, which 
is shown in the middle panel of Figure [1] fits the two 
bright peaks (MMS1 and MMS2) with an elliptical Gaus- 
sian and the third contiguous peak (MMS3) with a point 
source. After model subtraction, the residuals in the re- 
gion of the three main peaks are below 2 a, as shown in 
the right panel of Figure [TJ 

To compute the total masses of each of t hese contin- 
uum s ources, we follow the prescription of Hildcbrand 
((19831 ): 



where M, d, Si nt) k, T dust are the mass, distance, 
integrated flux, dust opacity per unit mass column 
density, and dust temperature of the source, respec- 
tively. B\ (Tdust) is the Planck function at wave- 
length A and temperat ure Tdust- We assume Tdust = 
30 K , d = 3.5 kpc (jHariu. Walmslev. fc Wouterloot] 
Il993h . a dust emissivity index of j3 = 1.5, k = 
0.1(250.0 iim/Xf g' 1 cm 2 (Hi ldebrandl [T983). and a gas- 
to-dust mass ratio of 100. Under these assumptions, we 
derive total masses of 225 M Q for MMS1, 205 M for 
MMS2, and 51 M fl for MMS3. 

IWilliams. Fuller, fc Sridharanl (|2004h observed 
IRAS 23033 using the JCMT and measured 850 /zm and 
450 /xm integrated fluxes of 10.4+0.2 Jy and 48.6+2.6 Jy, 
respectively, for a region encompassing MMS1, MMS2, 
and MMS3. Usin g these fluxes and assuming (3 = 0.9 
and Tdust = 52 K. (Williams. Fuller, fc Sridharanl (j2004D 
compute masses of 237 M at 850 fim and 106 M Q at 
450 /im for the IRAS 23033 complex. S imilarly, using 
single-d ish measurements at 1.2 mm, iBeuther et al.l 
(|2002al) derived a mass of 2327 M for IRAS 23033, 
assuming (3 = 2 and Tdust = 52 K. The differences be- 
tween the calculated masses are caused by three factors: 
the differences in the measured fluxes, the differences in 
the assumed dust opacities, and the differences in the 
assumed dust temperatures. Our interferometric maps 
filter out much of the flux visible in the single-dish maps. 
Because our measurements were taken at a different 
wavelength than any of the single-dish measurements, 
the computation of the exact amount of flux filtered out 
by the interferometer is dependent on the assumed dust 
opacity and temperature. 

If wc adopt the same flux-to-mass co nv ersion 
prescription as IWilliams. Fuller, fc Sridharanl (l200l 
(namely, assuming Td us t = 52 K and using the 
Osscnkopf fc Henningl (1994) dust opacities extrapolated 
to 3.3 mm) we calculate masses of 131 M , 120 M , 
and 29.5 M for MMS1, MMS2, and MMS3, respectively. 
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TABLE 1 
Molecular Lines Observed 



Frequency 


Molecular 


Array 


kA 


Velocity Res. 


rms / channel 


(GHz) 


Line 


Configurations 


Range 


(km s -1 ) 


(Jy beam -1 ) 


86.243440 


SiO 2-1 


B, C 


2-69 


0.337 


0.13 


86.340167 


H 13 CN 1-0 


B, C, D 


2-68 


0.085 


0.21 


86.754330 


H 13 CO+ 1-0 


B, C 


2-69 


0.337 


0.13 


89.188526 


HCO+ 1-0 


C, D 


2-25 


0.041 


0.75 


93.173777 


N 2 H+ 1-0 


B, C 


2-74 


0.079 


0.079 


96.412961 


C 34 S 2-1 


B, C 


2-74 


1.21 


0.058 


96.741377 


CH 3 OH 2-1 


B, C 


2-74 


0.302 


0.12 




23 h 05 m 27 s 26 s 25 s 24 s 23 s 22 s 26 s 25 s 24 s 23 s 22 s 26 s 25 s 24 s 23 s 22 s 

R. A. (J2000) R. A. (J2000) R. A. (J2000) 



Fig. 1. — Continuum emission at 3 mm (left), a model of the continuum emission consisting of 2 elliptical gaussians and a point source 
(middle), and the residuals obtained after subtracting the model from the data (right). The contours start at 3cr in the two left panels and 
2<r in the right panel; the interval between contours is la in all three panels, where a = 2.9 mjy beam -1 . The ellipse in the lower right of 
the left panel shows the size and orientation of the syn thesized beam. The star, triangle, and cross with error ellipse indicate the positions 
of the MSX source, H2O maser (Bcuthcr ct al. 2002b), and IRAS point source, respectively. 



We note, however, that the lOssenkopf &; Hennind (|1994f ) 
opacities were not originally tabulated for wavelengths as 
long as 3.3 mm, nor were they specifically meant to apply 
to massive sta r-forming regions. Sim ilarly, if we use the 
prescription of lBeuther et alj (|2002af) for converting inte- 
grated flux to mass, we find masses of 459 M , 419 M , 
and 103 M for MMS1, MMS2, and MMS3, respectively. 
Clearly the mass uncertainties are quite large. Neverthe- 
less, even allowing for generous uncertainties in the flux- 
to-mass conversion and in the star formation efficiency, 
at least MMS1 and MMS2 appear to have enough mass 
to each form one or more massive stars. 

Gravitationally bound cores capable of forming mas- 
sive stars typically have very high densities, n > 
10 5 cm -3 , and tend to have yirial r atios, ow = M v j r /M , 
close to unity (jPlume et al.lll997t IMcKee k. T an 2003). 
Do either of the spatially resolved cores in IRAS 23033 
meet these criteria? We can compute the mean densities 
of the two brightest cores using their total masses (as- 
sumed to be 225 M for MMS1 and 205 M for MMS2, 
as computed above) and their mean radii, defined as half 
of the mean of their major and minor axes, as deter- 
mined from the fits shown in Fig [TJ These mean radii 
are 0.14 ± 0.02 pc for MMS1 and 0.15 ± 0.03 pc for 
MMS2. In this way, we compute mean gas densities of 
(3 ± 1) x 10 6 cm -3 and (2 ± 1) x 10 6 cm -3 , for MMS1 



and MMS2, respectively. Such high gas densities again 
qualify both cores as potential sites of massive star for- 
mation. 

Figure [5] shows integrated intensity maps of 
IRAS 23033 in four molecular lines which princi- 
pally trace dense gas: N 2 H+ 1-0, H 13 CO+ 1-0, H 13 CN 
1-0, and C 34 S 2-1. These four dense gas tracers 
peak at different locations. The best column density 
tracers should be N2H" 1 " 1-0 and dust, because both of 
these should be optically thin and relatively (although 
perhaps not entirely) unaffected by the complex carbon 
chemistry which can change the abundances of H 13 CO + 
and H 13 CNThere is no indication in the N 2 H + loi - 0i2 
line shape that this transition is optically thick. Indeed, 
comparison of Figure [1] and Figure [2] shows that the 
spatial distribution of N-2H+ is very similar to the 
distribution of dust. The peak positions of N2H + and 
dust are separated by less than a beam width in both 
MMS1 and MMS2 (see also Fig. [3]). Hence, we conclude 
that, as expected, N2H + 1-0 and dust both trace the 
densest gas in the region. We can therefore use the 
N2H+ loi - 0i 2 line widths and the masses determined 
from the continuum emission to compute the virial 
ratios for MMS1 and MMS2. We derive the N 2 H+ line 
widths from the isolated (i.e. non-blended) hypcrfinc 
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component of the N2H + line. The N2H + loi-Oia line 
widths (full-width, half-power) in MMS1 and MMS2 
are 2.5 ± 0.3 and 2.4 ± 0.2 km s -1 , respectively. If we 
approximate the cores as spheres with uniform densities 
equal to the mean densities stated above and assume 
that the shape of the N2H + I01-O12 line is approximately 
Gaussian, then the virial masses of the cores are given 

by 



M vir = 



5R(Av) 2 
8Gln2 



(2) 



where R is the radius of the core and Av is the observed 
N 2 H + line width. By this prescription, the virial masses 
of MMS1 and MMS2 are 178 and 184 M respectively. 
Hence, their virial ratios are 0.79 and 0.90. The virial 
parameter can be interpreted as a measure of a cloud's 
degree of gravitational binding, with values less than 1 
indicating a gravitationally bound cloud. Therefore, sub- 
ject to the uncertainties in determining both the virial 
and dust continuum masses of the IRAS 23033 cores, we 
conclude that they are both gravitationally bound. This 
interpretation is consistent with the fact that MMS1 is 
clearly forming at least one star. It also suggests that 
MMS2 might already be forming a star, or may eventu- 
ally do so. 

Based on C 18 observations, lOnishi et al.l (|1996l ) 
derived a mean virial parameter of ~0.6 for star- 
forming molecular cloud cores in Taurus. Similarly, 
iTachihara. Mizuno. fe Fukuil (|2000D derive a mean virial 
ratio of 0.9 for C 18 cores in p Oph, with smaller ra- 
tios being associated with more actively star-forming 
cores. The C 18 cores in Taurus and p Oph have typ- 
ical masses in the tens of solar masses and radii of 
~0.2 pc. Despite the fact that these cores may have suf- 
ficient mass to form high-mass stars, the evidence from 
their environments is that they will instead form clus- 
ters of low-mass stars. The principal distinguishing fac- 
tor between the higher-mass cores in p Oph and those 
in IRAS 23033 seems to be that the IRAS 23033 cores 
have muc h higher surface densities: ~ 10 215 cm~ 2 
in p Oph (jTachihara. Mizuno. fc Fukuill200Ct ) vs. iV H . 2 ~ 
10 24 cm" 2 in IRAS 23033. The high surface densities of 
the IRAS 23033 cores are typical of those found in the 
massive star- f ormin g clumps surveyed in CS and C 34 S by 
iPlume et all (|1997f) . Thus, very dense gas seems to be 
key requirement for the formation of a massive star, as 
predicted by both the monolithic collapse and competi- 
tive accretion scenarios. 

3.1. Chemical Signatures of Evolutionary States 

The rapid evolution of the physical conditions in 
massive star-forming objects drives substantial chem- 
ical evolution in these objects. The time evolution 
of the abundances of particular molecules has not yet 
been well established for massive cores, but we can 
take some cues from the chemical studies of low-mass 
star-forming cores. In particular, we note that nu- 
merous studies have shown that, as a collapsing low- 
mass core heats up, CO desorption from grain man- 
tles causes substantial destruction of N2H" 1 ". The de- 
struction of N2H" 1 " lowers t he N2H" 1 " abundance near 
the center of the core (e.g. iLee. Bergin fc Evans 12004 



lDi Francesco. Andre, fc Mversll2004l ; iBergin et al.l 2002) 



This drop in the N2H + abundance has been observed 
in the centres o f several low-ma ss cores, such as L48 3 
(lJ0rgensenl2OO4D and Barnard lc (jMatthews et al.l2006h . 
iPirogov et al.l ( 2007f ) have shown that, in massive star- 
forming objects containing IRAS point sources, there is a 
measurable decrease in the N2H" 1 " abundance toward the 
column density peaks, again likely indicative of N2H + 
destruction in these evolved sources. 

Although the sequence of chemical processes at work in 
high-mass star formation likely differs from that in low- 
mass star formation, the destruction of N2H + by CO is 
likely to occur in both instances. In IRAS 23033, the dis- 
placement between the dust and N2H + emission peaks in 
MMS1, although small, suggests that the destruction of 
N2H" 1 " near the massive protostar in MMS1 may already 
be underway. Figure [3] shows our 3 mm continuum image 
overlaid with the peak positions of five of the molecular 
species we observed, as well a s the positions of the H2O 
mascr (taken from Table 2 of lBeuther et al.l l2002b). the 
MSX point source, and contours of 3.6 cm VLA radio 
continuum emission 1 . The coincidence of the radio con- 
tinuum and MSX emission with the dust peak in MMS1 
suggests that a massive protostar is forming there. The 
displacement between the dust and N 2 H + peaks suggests 
that heating by the massive protostar may have caused 
the desorption of CO from dust grain mantles and the 
subsequent destruction of N2H 4 " near the MMS1 peak. 
We emphasize, however, that the displacement between 
the N2H" 1 " and dust peaks in MMS1 is less than a beam 
width, and may not be significant. It is merely suggestive 
of chemical evolution in MMS1. 

Similarly, we can interpret the close correspondence of 
the dust and N2H" 1 " peaks in MMS2 as evidence that it 
is cold, chemically unevolved, and presumably therefore 
also very young. Its mass and apparent lack of chemical 
evolution suggests that MMS2 is a good candidate pre- 
stellar massive star-forming object. 

Our interpretation of MMS1 as chemically and dynam- 
ically more evolved than MMS2 is supported by the virial 
ratios of the two cores: lower in MMS1, indicating more 
gravitational binding, and higher in MMS2, indicating 
less binding and therefore perhaps an earlier evolution- 
ary state (see § [3|) . 



4. SPECTRAL LINE RESULTS: KINEMATICS 
4.1. Outflows 

Several previous studies have demonstra ted the pres- 
ence o f at least one outflow in IRAS 23033. iMao fc Zengl 
(2003) mapped the region in CO 3-2 using the 10 m 
Heinrich Hertz Telescope telescope and found over- 
lapping red- and blue-shifted outflow lobes emanating 
from a point not precisely coincident with the IRAS 
point source. Using BIMA array CO 2-1 observations, 
iBeuther. Schilke. fc Guethl (|2004l ) estimated the mass 
entrained in this outflow to be 119 M Q . Among the lines 
we observed, three are known outflow tracers: HCO + 
1-0, CH3OH 2-1, and SiO 2-1. These three lines give 
strikingly different pictures of the outflow activity in the 
IRAS 23033 region. 

1 The position of the M SX point source was taken from the 
MSX point source catalog (.Egan et a! J 2003) . The 3.6 cm radio 
continuum data were provided by H. Beuther. 
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Fig. 2. — Integrated intensity of four lines in IRAS 23033. Counter-clockwise from the top left panel, the lines are: N2H" 1 ", H 13 CO + , 
C 34 S, and H 13 CN. The contours start at 3<r in all panels and increase by steps of 4, 1, 1, and la, respectively. The N2H+ intensity was 
integrated over all seven hyperfine components of the line. The sizes and orientations of the four beams are indicated by the ellipses in the 
lower right corners of each panel. White crosses indicate the positions of the 3 mm continuum peaks from Fig. [l] The 1 <r rms in the panels 
are, from left to right, 290 mjy beam -1 km s — 1 , 210 mjy beam -1 km s -1 , 150 mjy beam -1 km s — 1 , and 320 mjy beam -1 km s — 1 . 



4.1.1. HCO+andSiO 

Figure 2] shows contours of red- and blue-shifted HC0 + 
from our BIMA array observations superimposed on the 
3 mm continuum map. This outflow clearly corresponds 
to the one seen by previou s auth ors (|Mao fc Zendl2003l : 
iBeuther. Schilke. fc Gueth|[200l . The HCO+ outflow 
is difficult to image because there is a significant back- 
ground of HC0 + associated with the underlying molec- 
ular cloud core, a source of confusion possibly amplified 
by the limited spatial resolution of our HCO + observa- 
tions. In constructing Figure^ we chose to display emis- 
sion originating only in those velocity ranges which best 
highlight the outflow emission. Curiously the SiO 2- 
1 emission in the region does not appear to trace the 
same outflow as seen in HCO + , or at least not the same 
part of it. Figure [5] shows the integrated intensity of 
SiO 2-1 over two velocity ranges: one centered at the 
systematic velocity and the other exactly matching that 
used to highlight the red lobe of the HCO + outflow in 
Figure |U We did not detect any SiO 2-1 emission at 
the 3cr level over the range of velocities which define the 
blue lobe of the HCO + outflow. The only suggestion of 
an association between the HCO + outflow and the SiO 
emission is that they are both extended and collinear. 



However, the velocity structures of the emission in the 
two molecules preclude their interpretation as parts of 
the same outflow. A more likely scenario is that the 
SiO peak at the systematic velocity is associated with 
hot and/or shocked gas in the immediate vicinity of the 
massive protostar. Figure [5] shows that the SiO emission 
peak at the systematic velocity is spatially coinciden t 
with the H2O maser identified by IBeuther et all (|2002bh . 
The two red-shifted SiO emission peaks to the northwest 
are harder to explain. These peaks may represent the 
red lobe of a separate outflow from the one delineated by 
the HCO + emission, or they may be entirely unrelated 
to outflow activity. If these SiO emission peaks are un- 
related to outflow emission-if, say, they were associated 
with other protostars-it is difficult to explain why they 
do not appear in other tracers. If the SiO peaks are asso- 
ciated with outflow activity, we expect that this outflow 
would be either weaker or chemically very different from 
the prominent outflow seen clearly in HCO + . 

The identity of the driving source of the HCO + out- 
flow is not clear. There is no candidate protostellar 
source at the point where the two outflow lobes overlap. 
The driving source is probably the MSX point source, 
which is spatially coincident with the peak of MMS1 and 
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Fig. 3. — Molecular line emission peaks (symbols) overlaid on 
the 3 mm continuum map (grey scale and contours). The mean- 
ing of each symbol is indicated in the legend. For simplicity, the 
positional error bars are not shown; they are typically ±1—2" in 
each coordinate. The thin white contours underlying the MSX 
point source a r e the V LA 3.6 cm radio continuum contours from 
IBeuther et"all [12002 pri , which have 0.7" resolution. The coinci- 
dence of the continuum and N2H+peaks in MMS2 suggests it is 
cold, and may be a pre-stellar massive star-forming core. The 
open ellipse in the lower left corner shows the size and orientation 
of the synthesized beam of the 3 mm continuum data. 
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Fig. 5. — Integrated intensity of SiO 2-1 plotted over 3 mm con- 
tinuum emission. The thick grey contours show the intensity of SiO 
2-1 integrated over the same velocity range as the red lobe of the 
HCO+ 1-0 outflow, namely -44.8 to -46.9 km s -1 (see Fig.gJ. 
The black contours show the intensity of SiO 2-1 around the sys- 



tematic velocity, integrated from —51 to —55 km s 



No SiO 



emission was detected over the velocity range corresponding to the 
blue lobe of the HCO" 1 " outflow. The thin grey contours and grey 
scale represent the continuum as in Fig. [4] All contours start at 3a 
and increase by steps of lcr, where a = 120 mjy beam -1 km s — 1 
for the red lobe and a = 150 mjy beam -1 km s -1 for the emis- 
sion integrated around the systematic velocity. The white triangle 
indicates the position of the water maser, which is coincident with 
the local peak of the SiO emission. The ellipse at the lower right 
shows the size and orientation of the SiO beam. 
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Fig. 4. — Integrated intensity of high-velocity HCO + 1-0 emis- 
sion plotted over continuum emission (grey scale and thin solid con- 
tours). The blue outflow lobe (dotted contours) is integrated over 
the velocity range —61.2 to —63.3 km s -1 and the red lobe (thick 
solid contours) is integrated over the range —44.8 to —46.9 km s — . 
Contours begin at 6 a for the red lobe and 4 <r for the blue lobe, 
where o = 108 mjy beam -1 km s — 1 , and increase by steps of 1 o. 
The ellipse in the lower right shows the size and orientation of the 
HCO+ beam. The blue-lobe HCO+ emission prominently overlaps 
MMS1, though their peaks do not coincide. The overlap of the 
red and blue lobes is probably an indication that the outflow is 
strongly inclined to the plane of the sky. 



the 3.6 cm continuum emission. This source is presum- 
ably the principal massive protostar in the region. Less 
clear, however, is whether the H2O maser and SiO emis- 
sion are also associated with this massive protostar, or 
whether they trace another source entirely. Just as the 
MSX source aligns neatly with the continuum peak of 
MMS1, the H2O maser aligns neatly with the peak of 
the systematic- velocity SiO emission. Neither pair of fea- 
tures coincides with the other. Hence, we speculate that 
there may be two massive protostars within the MMS1 
core, one situated at the millimeter continuum peak and 
one displaced a few arcseconds to the south. 

There is no unambiguous evidence of an outflow em- 
anating from MMS2. The HCO+ spectrum at the peak 
of MMS2 does show significant line wings (see Fig. |SJ), 
but these could be caused by 'spill-over' from the HCO + 
outflow emanating from the vicinity of MMS1. As can 
be seen in Figure [4] some of the HCO + emission which 
appears to be associated with the MMS1 outflow also 
overlaps MMS2. There is some evidence in the HCO + 
data cube of a second HCO + outflow in the region, but 
the spatial resolution of the data make it difficult to as- 
sociate it conclusively with either MMS1 or MMS2. To 
determine conclusively whether there is an outflow em- 
anating from MMS2, we would require data with higher 
spatial resolution and preferably in a transition which 
traces outflows more selectively than does HCO + . 

4.1.2. CH3OH 

As in many other high-mass star-forming regions, the 
CH3OH 2-1 emission in IRAS 23033 is very complex. 
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Fig. 6.— Integrated intensity of CH 3 OH 2-1 from -40.4 to - 
67.7 km s — 1 , plotted over 3 mm continuum emission. The black 
contours show the intensity of CH3OH 2—1 integrated over the two 
brightest components of the line, starting from 3<r and increasing 
by increments of 2tr, where a = 355 mjy beam -1 km s" 1 . The 
thin grey contours and greyscale represent the continuum as in 
Fig. [T] The size and orientation of the CH3OH synthesized beam 
are indicated by the ellipse in the lower right. The features labeled 
in the diagram are (1) a possible CH3OH counterpart to the SiO 
outflow, (2) a possible second outflow seen only in CH3OH, and 
(3) a conical structure which may represent the edges of a conical 
cavity carved out by a third outflow. Both of the outflows suggested 
by features (2) and (3) appear to lie approximately in the plane of 
the sky. 

Figure [6] shows the total integrated intensity of CH3OH 
in IRAS 23033. The CH 3 OH emission in IRAS 23033 is 
distributed widely in velocity and, in general, does not 
correlate neatly in either velocity or position with other 
outflow features. The exception is the elongated feature 
labeled '1' in Figure which is coincident both spatially 
and in velocity with the SiO peaks shown in Figure \5[ 
A prominent bar of CH3OH emission overlaps most, but 
not all, of the continuum-emitting region. The emission 
peak of the CH3OH bar does not coincide with either 
continuum peak. 

Although the CH 3 OH emission does not correlate well 
with the HCO+ outflow, it can be interpreted as tracing 
several other possible outflows in the region. For exam- 
ple, the dashed line labeled '2' in Figure [6] connects the 
CH3OH peak at the eastern edge of the image through a 
ridge of CH3OH emission to a second CH3OH peak near 
the western edge of the image. In velocity, the CH3OH 
peaks at the two extremes of line '2' are separated by only 
a few km s _1 but their alignment suggests they may be 
part of an outflow lying nearly in the plane of the sky. 
The lack of detected continuum emission at the position 
of either CH3OH peak supports this interpretation. The 
driving source of this possible outflow is not evident. 

The most prominent CH3OH emission feature is the 
large parabolic structure, labeled '3', which dominates 
the southern half of the image. This feature has no 
counterpart in any other tracer we observed. Its ve- 
locity structure suggests it may be the conical cavity 
of an outflow lying approximately in the plane of the 
sky. Again, the source of this possible outflow is not evi- 
dent, although the alignment of the feature is consistent 



with a driving source within the MMS2 core. If the driv- 
ing source did lie within MMS2, the second outflow lobe 
could be confused with the elongated bar of emission at 
the systematic velocity traced by the continuum emission 
and in several of the lines. 

4.2. A "Flattened Rotating Object" 

Four of the molecular transitions we observed have 
high critical densities and are known to trace dense 
star- forming gas. These are the transitions for which 
integrated intensity maps are shown in Figure O All 
four transitions trace similar structures in IRAS 23033: 
roughly elliptical distributions with long axes running 
from northwest to southeast, encompassing regions of 
high column density as traced by continuum emission 
(see Fig. [TJ The long axis of this elongated structure 
runs perpendicular to the axis of the prominent HCO + 
outflow. A fit to the H 13 CO + integrated intensity map 
shows that the position angle of the elongated feature 
is 35°±5° degrees, while that of the HCO + outflow is 
approximately 130°, meaning that the position angles of 
the two features differ by about 95°. 

As shown in Figure El the four dense gas tracers all 
peak at slightly different locations. As we discussed in 
§ [31 N2H + and dust both trace the regions of highest 
column density. Our C 34 S detection is too weak to de- 
termine the peak position accurately, except to say that 
it appears to correspond to the peak of the MMS1 core. 
The peak of the H 13 CN emission matches that of the SiO 
emission and the position of the H2O maser which, as we 
described in § 14.1.11 is displaced by several arcseconds 
from the coincident dust continuum, radio continuum, 
and MSX peaks. Notably, however, the maser and the 
emission peaks of H 13 CN SiO, and C 34 S all lie within 
one beam width of the MMS1 continuum peak, implying 
an association with MMS1. 

The remaining dense gas tracer, H 13 CO + , shows a 
distribution that requires some explanation. H 13 CO + 
peaks strongly at a point equidistant between the MMS1 
and MMS2 continuum peaks, in a region of substantially 
lower column density as indicated by the dust emission. 
The position of the H 13 CO + peak does not match that 
of any other feature in the region (see Fig. [3]). It is 
unlikely that this peak is associated with any coherent 
physical structure; it is probably an artifact of deple- 
tion of H 13 CO+ onto dust grains at the column density 
peaks, leaving a remnant "peak" where depletion effects 
are least pronounced. Although N2H + 1-0 and H 13 CO + 
1-0 peak at different locations, their kinematics clearly 
reveal that they trace the same contiguous structure in 
which the MMS1 and MMS2 peaks are embedded. The 
axis ratio of the elongated structure seen in H 13 CO + is 
3.1. Figure [7] shows position- velocity slices in N 2 H + Iqi- 
O12 and H 13 CO + taken along the major axis of the elon- 
gated structure, at a position angle of 35°. As shown 
in the right panel of Figure [71 both transitions exhibit 
the same velocity gradient, meaning that they trace the 
same large-scale structure. Note that the signal-to-noise 
ratio of the H 13 CN and C 34 S data is too low to permit a 
meaningful analysis of their velocity structures. 

On large spatial scales, the elongated structure seen 
in Figure [2] shares many characteristics in common with 
a stellar accretion disk. For example, its long axis lies 
perpendicular to at least one outflow, there is a veloc- 




velocity (km/s) velocity (km/s) velocity (km/s) 

Fig. 7. — Position-velocity slices taken along the long axis of the elongated structure seen in N2H+ I01-O12 (left) and H 13 CO + 1-0 
(center). The position angle of the position-velocity cut is 35°. Contours begin at 3<r and increase by intervals of la for N2H+ lot— O12 and 
2(7 for H 13 CO+. The white triangl es in the left and center panels show the position-centroided velocities of the emission in each channel. 
The right panel plots the two position-velocity curves on the same scale, illustrating the similarity between them. 



ity gradient along its major axis, and it appears only in 
the continuum and emission from molecules which trace 
dense, non-outflowing gas. However, the structure also 
has several characteristics which preclude its interpre- 
tation as an accretion disk, the first of which is its size. 
The major axis of the H 13 CO + emission over which there 
is a clear velocity gradient spans approximately 25", or 
about 87,000 AU at 3.5 kpc. By comparison, proto- 
stellar accretion disks around low- mass stars typically 
have radii in the ran ge 50-500 AU (|Dutrev et all 119961 ; 
iKitamura et al.ll2002f ). Little is known about the plausi- 
ble range of radii for accretion disks around massive pro- 
tostars. Table [2] gives an (incomplete) summary of the 
reported properties of candidate disks around massive 
protostars. In most cases, the evidence that these objects 
are disks consists of three points: they are elongated in 
either the continuum or a dense gas tracer, the axis of 
elongation is approximately perpendicular to that of a 
detected outflow, and there is a velocity gradient along 
the long axis of the object observed in either molecular 
line or maser emission. We note, however, that the range 
of reported physical properties of these disks is very large. 
For example, their radii range over more than two orders 
of magnitude. If we were to submit IRAS 23033 as a can- 
didate for this list based on its outflow properties and the 
shape and kinematics of its dense gas feature, it would 
only top out the list of an already very diverse sample. 
We note also that, while the molecular gas distribution in 
H 13 CO + and possibly N2H + could support the disk inter- 
pretation, the dust distribution does not. The dust emis- 
sion is better described as three strong peaks whereas an 
accretion disk should produce a centrally peaked distri- 
bution unless it had substantially fragmented. 

How, then, should we interpret the kinematics of 
IRAS 23033 and, more generally, of very large struc- 
tures with rotation signatures surrounding massive pro- 
tostars? We suggest that the observations to date, in- 
cluding those summarized in Tabled are consistent with 
the hypothesis that massive protostars may be fed by Ke- 
plerian accretion disks (e.g. IRAS 20126+4104), but that 



these accretion disks are smoothly connected to much 
larger, flattened, rotating structures. These "flattened 
rotating objects" are likely simply the remnants of the 
molecular cloud fragments from which the massive stars 
formed. Existing observations, typically made with sin- 
gle interferometer configurations, only reveal rotational 
motions on the spatial scales visible with that interfer- 
ometer configuration. Confirmation of the hypothesis 
that these rotating structures may extend from tens to 
thousands of AU would require multiple interferomet- 
ric observations of each source with different interferom- 
eter configurations. Such observations exist in only a 
few cases. IRAS 201 26+4104 is an instructive example: 
ICesaroni et al.l (|1999l ) detect apparently Keplerian rota- 
tion on scales down to a few 100 AU using high-resolution 
CH3C N observations while IZhang. H unter. &: Sridharanl 
(|1998fl detect rotation out to radii of 5,000 AU in NH 3 . 
We interpret our observations of IRAS 23033 in a sim- 
ilar manner. The optically thin H 13 CO + line depicts 
a contiguous, flattened, slowly-rotating structure with a 
radius as large as 40,000 AU. N2H 4 " and dust continuum 
emission respectively reveal two and three condensations 
of higher column density within this apparently contigu- 
ous structure. The cores detected in N 2 H + and dust 
continuum emission may be fragmentation products of 
the original molecular cloud from which the IRAS 23033 
complex formed. 

It is also possible that the velocity gradient seen in 
H 13 CO + and N 2 H + may simply indicate that the elon- 
gated structure is a sheared filament. This hypothesis 
cannot be readily ruled out using the data in hand. 

4.3. Infall 

iZhoul (|1992D showed that, in a dense molecular 
cloud core undergoing some form of inside-out col- 
lapse, optically thick lines will be skewed to the blue 
by an amount which depends on the infall velocity. 
The combination of a double-peaked, asymmetrically 
blue HCO + line whose self-absorption peak is cen- 
tered at the velocity of the optically thin H 13 CO + 
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TABLE 2 

Candidate Massive Protostellar Disks 



Object 


Distance 


Radius 


Disk Mass 


Ref. 


(kpc) 


(AU) 


(M Q ) 





G192. 16-3.82 2 130 

Cepheus A HW2 0.73 330 
Becklin-Neugebauer Object 0.46 800 

AFGL 5142 1.8 900 

IRAS 18089-1732 3.6 ~1000 

NGC 7538 IRS IN 2.8 1000 

IRAS 20126+4104 1.7 1700 

G24.78 Al 7.7 4000 

G24.78 A2 7.7 4000 

G24.78 C 7.7 8000 

G31.41 7.9 8000 

NGC 7538 S 2.8 15,000 



3 
1-8 



~16 

10 
23 
23 
23 
23 
400 



Kurtz (2001 



Shepherd. Clausscn, 
Patel et al. (2005) 
Jiang (2005) 
Zha ng et al. (2002) 
Beuther et al. (2005') 
Pestalozzi"eTal^2ro4 ) 
Cesaroni et al. (,1997, 1999) 
Beltran et ai"7200lT " 
Belt ran et al. (2004) 
Belt ran et al. (2004) 
Beltran et al. (2004) 
Sandell. Wright, fc Forster (2003) 



line has been interpreted as an indica tor of infall in 
both low-mass dGregersen fc Evans!l2000l) and high -mass 
((Fuller. Williams, fc Sridharanll2005D~ star-forming cores. 
Figure M shows spectra of HCO+ and H 13 CO+ at the 
three continuum peaks seen in Figure [TJ At the contin- 
uum peaks of both MMS1 and MMS2, the HCO+ line ap- 
pears self-absorbed, indicating that it is optically thick. 
At the position of MMS1, the line is slightly skewed to 
the red. We note, however, that there is significant out- 
flow contamination in the spectrum of the core at this 
position. At the position of MMS2, the line profile is 
skewed to the blue, suggesting it may be collapsing. 

Another possible interpretation of the double-peaked 
HCO + line profiles is that there are two sources along 
the line of sight. This possibility would seem to be ruled 
out by the fact that the optically thin H 13 CO + line is 
single-peaked throughout the region and it always peaks 
at the velocity of the HCO + self-absorption feature. A 
third possible interpretation for the double-peaked and 
asymmetrically blue line profiles is that emission from 
the two cores is being blurred together in the compara- 
tively large (~ 10") HCO + beam. To assess the likeli- 
hood of this possibility, we smoothed the H 13 CO + data 
cube so that it had the same beam orientation and size 
as the HCO + beam. Even in the smoothed data set, the 
H 13 CO + line remains single-peaked at each continuum 
peak, thereby eliminating resolution effects as a likely 
explanation for the asymmetric line profiles. Hence, we 
conclude that there is evidence for collapse in MMS2. 

5. SUMMARY 

We used millimeter intcrfcromctry to study the struc- 
ture and kinematics of the massive star-forming region 
IRAS 23033+5951. We find that the region resolves into 
at least three distinct cores, at least two of which have 
sufficient mass to form one or more massive stars. There 
is strong evidence from the spatial coincidence of strong 
radio continuum emission, an MSX point source, and a 
3 mm continuum peak that MMS1 is forming a single 
massive star. This core also powers at least one outflow 
seen in HCO + 1-0 emission and possibly a second outflow 
seen in SiO 2-1 and CH3OH 2-1. The low virial ratio 
and moderately high density of MMS1 suggest it could 
be forming a massive star by collapse and disk accre- 
tion according to the "monolithic collapse" model. The 
distribution of N2H + in IRAS 23033 strongly resembles 
that of the dust emission, suggesting that N2H+ is trac- 
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Fig. 8.— Spectra of HCO+ {thick black lines) and H 13 CO+ 
(thin black lines) at the positions of MMS1 (top), MMS2 peak 
(middle), and MMS3 (bottom). The vertical line in each plot in- 
dicates the fitted peak of the H 13 CO + spectrum, an indication of 
the systematic velocity at that position. 



ing the locations of highest column density. Conversely, 
the fact that H 13 CO + peaks between the column density 
peaks revealed by the dust and N2H + emission probably 
reflects the depletion of H 13 CO + onto grain mantles in 
the dense regions. 

MMS2 appears to be quiescent: it harbors no maser or 
MSX sources and has no unambiguous outflow signature. 
The spatial coincidence of its N2H + and dust continuum 
emission peaks suggest that it is a cold core in an early 
stage of evolution. With an estimated total mass (dust 
and gas) of 205 M Q , a virial ratio of 0.9, and a possible 
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spectral signature of collapse, MMS2 appears poised to 
form a massive star of its own. It may be an example of 
a "Class 0" massive protostar. 

Both of the cores described above are embedded in 
a smooth structure seen most clearly in the optically 
thin H 13 CO+ 1-0 line, but also seen in H 13 CN 1-0 and 
HCO + 1-0. This structure has an axis ratio of 3.1, and 
its long axis lies almost exactly perpendicular to that of 
the HCO+ outflow. The H 13 CO+ emission peak of the 
structure is well-centered between the continuum peaks 
of the two bright cores. There is a clear velocity gradi- 
ent along the structure's long axis which is suggestive of 
rotation, although the size of the structure precludes its 
interpretation as a classical protostellar accretion disk. 
We suggest it may either be the rotating, flattened rem- 
nant of the natal molecular cloud from which the two 



Beltran, M. T., et al. 2004, ApJ, 601, L187 

Bergin, E. A., Alves, J., Huard, T., & Lada, C. J. 2002, ApJ, 570, 
L101 

Beuther, H., Schilke, P., Menten, K. M., Motte, F., Sridharan, 

T. K., & Wyrowski, F., 2002a, ApJ, 566, 945 
Beuther, H., Walsh, A., Schilke, P., Sridharan, T. K., Menten, 

K. M., & Wyrowski, F. 2002b, A&A, 390, 289 
Beuther, H., Schilke, P., Sridharan, T. K., Menten, K. M., 

Walmsley, C. M., & Wyrowski, F. 2002c, A&A, 383, 892 
Beuther, H., Schilke, P., & Gueth, F. 2004, ApJ, 608, 330 
Beuther, H., Zhang, Q., Sridharan, T. K., & Chen, Y. 2005, ApJ, 

628, 800 

Bonnell, I. A., Bate, M. R., & Zinnecker, H. 1998, MNRAS, 298, 
93 

Bonnell, I. A., Bate, M. R., Clarke, C. J., & Pringle, J. E. 2001, 

MNRAS, 323, 785 
Bonnell, I. A., Clarke, C. J., Bate, M. R., & Pringle, J. E. 2001, 

MNRAS, 324, 573 
Bonnell, I. A., Vine, S. C, & Bate, M. R. 2004, MNRAS, 349, 735 
Cesaroni, R., Felli, M., Testi, L., Walmsley, C. M., & Olmi, L. 1997, 

A&A, 325, 725 

Cesaroni, R., Felli, M., Jenness, T., Neri, R., Olmi, L., Robberto, 

M., Testi, L., & Walmsley, C. M. 1999, A&A, 345, 949 
Di Francesco, J., Andre, P., & Myers, P. C. 2004, ApJ, 617, 425 
Dutrey, A. et al. 1996, A&A, 309, 493 

Egan, M. P. et al., 2003, in Air Force Research Laboratory 

Technical Report AFRL-VS-TR-2003-1589 
Williams, S. J., Fuller, G. A., & Sridharan, T. K. 2004, A&A, 417, 

115 

Fuller, G. A., Williams, S. J., & Sridharan, T. K. 2005, A&A, 442, 
949 

Gregersen, E. M. & Evans, N. J. 2000, ApJ, 538, 260 
Harju, J., Walmsley, C. M., & Wouterloot, J. G. A. 1993, A&AS, 
98, 51 

Hildebrand, R. H. 1983, QJRAS, 24, 267 
Jiang, Z., et al. 2005, Nature, 437, 112 
j0rgensen, J. K. 2004, A&A, 424, 589 



current cores formed, or simply a sheared filament. 



M. A. Reid is currently supported by a fellowship from 
the Harvard-Smithsonian Center for Astrophysics but, at 
the time the data for this paper were acquired, he was 
supported by a scholarship from the Natural Sciences 
and Engineering Research Council of Canada. B. C. 
Matthews is supported by the National Research Coun- 
cil of Canada. The BIMA array was operated with sup- 
port from the National Science Foundation under grants 
AST 02-28963 to UC Berkeley, AST 02-28953 to U. Illi- 
nois, and AST 02-28974 to U. Maryland. FCRAO is 
supported by NSF grant AST 02-28993. The radio con- 
tinuum data shown in Figure [3] were generously provided 
by H. Beuther. 



Kitamura, Y. et al. 2002, ApJ, 581, 357 

Krumholz, M. R., McKee, C. F., & Klein, R. I. 2005, Nature, 438, 
332 

Krumholz, M. R. 2006, ApJ, 641, L45 

Lee, J.-E., Bergin, E. A., & Evans, N. J. 2004, ApJ, 617, 360 
Mao, R. Q., & Zeng, Q. 2003, AcASn, 44, 61 

Matthews, B. C, Hogerheijde, M. R., j0rgensen, J. K, & Bergin, 

E. A. 2006, ApJ, 652, 1374 
McKee, C. F., & Ostriker, E. C. 2007, ARA&A, 45, 565 
McKee, C. F. & Tan, J. C. 2003, ApJ, 585, 850 
Onishi, T., Mizuno, A., Kawamura, A., Ogawa, H., & Fukui, Y. 

1996, ApJ, 465, 815 
Ossenkopf, V. & Henning, Th. 1994, A&A, 291, 943 
Patel, N., et al. 2005, Nature, 437, 109 

Pestalozzi, M. R., Elitzur, M., Conway, J. E., & Booth, R. S. 2004, 
ApJ, 603, L113 

Pirogov, L., Zinchenko, I., Caselli, P., & Johansson, L. E. B. 2007, 
A&A, 461, 523 

Plume, R., JafTe, D. T., Evans, N. J., Martm-Pintado, J., & Gomez- 
Gonzalez 1997, ApJ, 476, 730 
Sandell, G., Wright, M., & Forster, J. R. 2003, ApJ, 590, L45 
Sault, R. J., Teuben, P. J., & Wright, M. C. H. 1995, in ASP Conf. 
Ser. 77, Astronomical Data Analysis Software and Systems IV, 
ed. R. A. Shaw, H. E. Payne, & J. J. E. Hayes (San Francisco: 
ASP), 433 

Shepherd, D. S., Claussen, M. J., & Kurtz, S. E. 2001, Science, 
292, 1513 

Sridharan, T. K., Beuther, H., Schilke, P., Menten, K. M., & 

Wyrowski, F. 2002, ApJ, 566, 931 
Tachihara, K., Mizuno, A., & Fukui, Y. 2000, ApJ, 528, 817 
Zhang, Q., Hunter, T. R., & Sridharan, T. K. 1998, ApJ, 505, L151 
Zhang, Q., Hunter, T. R., Sridharan, T. K, & Ho, P. T. P. 2002, 

ApJ, 566, 982 
Zhou, S. 1992, ApJ, 394, 204 



